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ABSTRACT: We propose a definition for the entanglement entropy of a gauge theory on a 
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both Abelian and Non-Abelian gauge theories. For Zjv and U(l) theories, without mat¬ 
ter, our definition agrees with a particular case of the definition given by Casini, Huerta 
and Rosabal. We also argue that in general, both for Abelian and Non-Abelian theories, 
our definition agrees with the entanglement entropy calculated using a definition of the 
replica trick. Our definition, however, does not agree with some standard ways to mea¬ 
sure entanglement, like the number of Bell pairs which can be produced by entanglement 
distillation. 
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1 Introduction 

Entanglement entropy is an important measure of quantum correlations in a system and 
provides a way to quantify the “essential weirdness” which arises in quantum mechanics. 
Gauge theories are a central paradigm in modern physics. Three of the four fundamental 
forces in nature are described by gauge theories. And via holography, they are important 
now in our study of the fourth fundamental force, gravity. Also, gauge theories are an 
increasingly important corner stone of condensed matter physics. It is therefore important 
to try and provide a precise definition of entanglement entropy for gauge theories. 

This is however not straightforward. In a theory, like a spin system or a scalar held 
theory, which has localised physical excitations, the Hilbert space of states admits a tensor 
product decomposition in terms of local degrees of freedom. To obtain the entanglement 
entropy associated with a region of space we first trace over the degrees of freedom outside 
to obtain a density matrix. And then calculate the von Neumann entropy of this density 
matrix to get the entanglement entropy. In a gauge theory, unlike the spin system or the 
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scalar field theory, however, the physical degrees of freedom are not always local in space- 
time. For example, in a pure gauge theory without matter, they correspond to closed loops 
of electric or magnetic flux, which are non-local. As a result, the Hilbert space of physical 
gauge-invariant states does not admit a local tensor product description, and without such 
a decomposition it is not clear how to define a density matrix associated with some region 
of space and obtain the entanglement entropy. 

In this note we attempt to give such a definition for gauge theories. Since the entangle¬ 
ment entropy is known to suffer from short-distance divergences in field theories, we work 
with gauge theories living on a spatial lattice. Motivated by the work of Casini, Huerta 
and Rosabal, [1], (CHR), see also [2], we then ask a slightly different question from the one 
posed above. Instead of a spatial region we consider a subset of links in this spatial lattice 
and ask about the entanglement entropy of these links. We find that this question lends 
itself to a precise and elegant answer. 

The central point leading to our definition is as follows. While the Hilbert space 
of gauge-invariant states does not admit a tensor product decomposition in terms of local 
degrees of freedom, a larger Hilbert space T~L can be defined admitting such a decomposition. 
This larger Hilbert space is obtained by taking the tensor product of the Hilbert spaces on 
which the individual link excitations live, and therefore arises quite naturally in the gauge 
theory. It is larger because it includes both gauge-invariant and non-gauge-invariant states. 
By regarding any gauge-invariant pure state |i p) G 7i we can then obtain a density matrix, 
Pin, for the links of interest, which we call the “inside” links, by tracing over the outside 
ones. Our definition of the entanglement entropy is then simply given by taking the von 
Neumann entropy of this density matrix: 

See = ~ Tr Hin p in log p in , (1.1) 

where the trace is taken over the Hilbert space, 'Hi n , which is the tensor product of the 
Hilbert spaces of the inside links. 

The definition above applies to both Abelian and Non-Abelian gauge theories. Also, 
while for simplicity we will focus on 2 spatial dimensions in this paper, the definition holds 
in higher spatial dimensions as well. 

We find that in the Abelian Zjy and 1/(1) gauge theories, without matter, our definition 
agrees with the electric centre choice for the entanglement entropy given in [1], For the Non- 
Abelian case we find things are similar, and our definition can also be related to an electric 
centre type prescription, although the different sectors do not correspond to different values 
of electric flux. 

Importantly, we also find that for both the Abelian and Non-Abelian cases our definition 
agrees with the result for the entanglement entropy obtained by using a definition of he 
replica trick. 

Our definition has one unsatisfactory feature, though. In quantum information theory 
the entanglement of a bipartite system, i.e. a system with two parts, is often quantified 
by comparing it to the entanglement of a standard system, e.g., a collection of Bell pairs. 
The comparison can be done, by converting the entanglement between the two parts into 
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the entanglement of a set of Bell pairs, as happens in entanglement distillation. Or by 
consuming the entanglement in a set of Bell pairs to quantum teleport one of the two parts, 
as happens in entanglement dilution. It turns out, that the entanglement entropy we have 
defined does not equal the number of Bell pairs which can be produced or consumed in this 
way. This fact was noticed by CHR already for the Abelian case, and it continues to be 
also true, with an interesting new aspect, in the Non-Abelian case. 

The underlying reason for this difference is the presence of non-local gauge-invariant 
operators, and excitations associated with these operators in gauge theories. The protocols 
for entanglement distillation or dilution involve only local operations (and classical com¬ 
munication) and this turns out to be not enough for converting the full entropy into Bell 
pairs. Since the presence of non-local gauge-invariant operators is a central feature in gauge 
theories, as was emphasised at the outset above, we suspect that such a limitation might be 
more generally true in other attempts to define the entanglement entropy in gauge theories 
as well. 

The paper is structured as follows. We start first by discussing the Z 2 gauge theory 
without matter in 2 spatial dimensions in section 2. Some general properties for our defini¬ 
tion of the entanglement entropy are discussed in the context of the Z 2 example in section 
2.1. Then we consider the 17(1) case in section 3. This sets the stage for the Non-Abelian 
case, which we discuss specifically for the SU ( 2 ) example in section 4. Connections to the 
replica trick are analysed in section 5, followed by a discussion of additional aspects includ¬ 
ing the relation to more operational definitions of entanglement, in section 6 . We end in 7 
with some conclusions and open directions. 

Before ending the introduction let us also mention some of the other relevant literature. 
Background material on entanglement in held theories that we found useful is in [3], [4], [5], 
[ 6 ], [7], [ 8 ], and [9]. For issues related to entanglement in gauge theories some references 
include [10], [11], [12], [13], [14], [15], [16], [17], 

2 The Z 2 Lattice Gauge Theory 

Our starting point is a Z 2 lattice gauge theory without matter, one of the simplest examples 
of a gauge theory. For definiteness we work in 2 + 1 dimensions. We will consider the theory 
in the Hamiltonian formulation, where time is continuous, and the two spatial directions 
are on a lattice, see [19], which for simplicity we take to be square. 

A state in this theory is defined, at an instant of time, on the two-dimensional spatial 
lattice. There is a two dimensional Hilbert space, corresponding to one qubit, associated 
with each link of this lattice. A basis for this Hilbert space can be taken to be | ± 1), with 

<73 1 ± 1 ) = ± 1 | ± 1 ) ( 2 . 1 ) 

where <73 is the standard Pauli operator. This basis is orthonormal, meeting conditions 
(+1| + 1) = 1, (+1| — 1) = 0, etc. Gauge transformations act on vertices of the lattice. A 
gauge transformation Qy t acting on vertex V) flips the values of all link variables emanating 
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Figure 1. The subset of links represented by solid lines constitute the inside. The remaining links 
shown by dashed lines are the outside. 


from Vi. In the basis introduced above we have 

G Vi = n*i ( 2 - 2 ) 

where the product is taken over all links which emanate from V{. Gauge-invariant states 
are invariant under all gauge transformations, i.e., if |V’) is a gauge-invariant state, then, 

QVi\^) = \^) (2.3) 

for all gauge transformations Q\r. Similarly gauge-invariant operators are invariant under 
gauge transformations. If O is a gauge-invariant operator, 

QvfiQvl = 6. (2.4) 

Gauge-invariant operators in the theory include Wilson loops which are the products of 
<73 operators for links forming closed plaquettes, and the operators on each link. These 
measure the discrete analogue of the magnetic and electric fluxes respectively in this theory. 
The most general gauge-invariant operator can be constructed by taking products of these 
operators. In the discussion below we refer to the a\ operators, somewhat loosely, as the 
electric flux operators. 

Now consider a gauge-invariant state, | ip), and some subset of all the spatial links. We 
will call these links as lying in the “inside”, although our discussion is general and does not 
require the collection of links to form a closed or even connected region. For a concrete 
example, consider Fig 1 with the subset of links of interest lying inside and on the boundary 
of the rectangle. 

As was discussed in the introduction there is an extended Hilbert space, Ti, which is 
obtained in this case by taking the tensor product of all the two dimensional Hilbert spaces 
defined at each link. This Hilbert space by definition admits a tensor product decomposition 
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in terms of the Hilbert spaces of links lying inside, 7~Li n , and outside (i.e. the rest of the 
links), 7-Lout • A basis for 7~L is given by taking the tensor product of the |±1) basis introduced 
above for each link. Since the set of gauge-invariant states is a subset of 7~L, any state \ift) is 
automatically embedded in 7~L. It is easy to see that the corresponding state in 7~L is unique. 
For example, there is a unique way to expand \tjj) in terms of the basis obtained by taking 
the tensor product of the | ± 1) states on each link. We will sometime refer to this as the 
natural embedding of | ip) E 7~L below. 

Our proposal is to take the state \ip), now regarded as an element of 7~L, and define the 
density matrix pi n for the inside set of links by taking the trace, 

Pin = ( 2 - 5 ) 

The entanglement entropy is then given by the von Neumann entropy associated with pi n , 
eq.(l.l). It is clear that this definition can be applied to any set of links. 

Before proceeding let us also comment on some of the existing literature here. A similar 
definition for the Abelian case was given in [13]. Also, in the context of an extended lattice 
construction, as we will discuss further in section 6.3, the same definition was given, both 
for the Abelian and Non-Abelian cases, in 1 [12], [14], 

2.1 Some Properties 

Some comments are now in order about the properties for the entanglement entropy which 
follow from our definition. Even though we state them in the Z 2 context here, these 
properties are actually generally true for all Abelian and Non-Abelian cases we study in 
this paper. 

First, it is worth emphasising that the definition is unambiguous. Any gauge-invariant 
state corresponds to a unique state in 7~L, as was mentioned above, and for any state in 7i 
there is a unique density matrix pi n , leading to a unique entanglement entropy. In fact, 
one can similarly define the Renyi entropies for any set of links and these too are all then 
unambiguously defined. 

Second, we note that the density matrix defined above pi n allows one to calculate the 
expectation value of any gauge-invariant operator which only acts on the inside links. Any 
such operator, O, can be promoted to be an operator acting on 7~L. And if it only acts on 
the inside links then it is easy to see that 

{ip\0\4’) = Tr (pinO). (2.6) 

In this way we see that pi n contains all the information needed to describe the results of any 
physical measurement done on the inside links. Later on, in section 6.1, we will see that the 
converse is also true, in the following sense. Namely, we will argue that measurements done 
on the inside links (more correctly an ensemble of inside link systems) can fully determine 
Pin■ 

1 We thank W. Donnelly for correspondence in this regard. 
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Third, the Hilbert space B obtained from the tensor product of the two-dimensional 
Hilbert spaces at each link also inherits from these 2-dimensional subspaces a natural inner 
product, meeting all the standard properties of positivity etc. Thus the entanglement 
entropy we have defined automatically meets the strong subadditivity condition. If A, B 
and C are three sets of links which share no links in common then we have that 

Saub + Scub > Sb + SaubuC (2-7) 

where Saub is the entanglement entropy for the set of links, AU B etc. 

While these three properties have been discussed in the context of the Z 2 theory in this 
section, they will also continue to to be true for all the other theories considered in this 
paper essentially for the same reasons. 


2.2 Equivalence With Other Definitions 


Here we show that the definition for entanglement entropy given above by embedding the 
state in the Hilbert space B agrees with the electric centre choice of [1] (see also [2] for a 
clear discussion). 

To establish the equivalence we will need to re-express the entanglement entropy, 
eq.(l.l), as a sum over different, suitably defined, sectors. First, let us introduce some 
terminology. Vertices in the lattice will be said to lie inside if all links ending on them are 
inside links, and outside if all terminating links are outside links. In addition there will be 
boundary vertices which have some terminating links lying in the inside subset and some 
lying in the outside subset. 

To carry out the trace over Bout and define pi n we can work in any basis. For connecting 
to the electric centre prescription of CHR it is useful to consider the basis in which the 
outside links ending on a boundary vertex are in eigenstates of the electric flux operators 
defined on these links. For example, in Fig 2 we have shown the boundary vertex V\ which 
has two inside links, labelled as 3,4, and two outside links, labelled as 1, 2, ending on it. In 
the basis we work with, the two outside links are in eigenstates of the o\ operators. This is 
of course not a complete specification of the basis of B ou t • But the remaining steps below 
will not depend on specifying the basis for the other links in B ou t explicitly, so we do not 
do so here. 

Now at any boundary vertex let us define the total electric flux operator going outside to 
be the product of all the o\ operators acting on the outside links ending at this vertex. For 
example in Fig 2 vertex V\ has two outside links, labelled 1, 2, with electric flux operators 
which we denote to be respectively. The total electric flux operator going outside 

at this vertex is then 


a 


(!)_( 2 ) 


1 


<7 


1 


( 2 . 8 ) 


The basis we are working with therefore consists of eigenstates of the total out-going electric 
flux operator at each vertex. We can then define a vector, k, whose entries, k\, & 2 , ■ ■ ■ , give 
the eigenvalues of the total electric flux operators going outside at the boundary vertices 
Vi, V 2 , ■ ■ ■, respectively. This vector specifies the electric flux which leaves the inside region 
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Figure 2. The lattice with boundary vertex Vj . Links 1 and 2 are in the outside subset, while 
links 3 and 4 are in the inside subset. 


at all the boundary vertices. We note that for the Z 2 theory each entry in k takes values 

± 1 . 

It is now clear that the Hilbert space of the outside links, 'Houty can be written as a 
direct sum over sectors of definite electric flux, 



'Hout - 


(2.9) 

By expressing the identity operator, 1, of 'H.outy as a sum over projection operators P k into 
the subspace 

1 = E i3k , ( 2 - 10 ) 

k 

we can express pi n as 

Pin = TrHoJ^K^I =^2^Hout 

k 

(p k |V.)(V-l) ■ 

(2.11) 

So that 

Pin — Pin ? 

k 


(2.12) 

with, 

pf n = Tr W kJV’)(</i 

(2.13) 

where, 

= TrHout ( P k 

IV’XV’l) • 

(2.14) 


Any outside vertex by definition has only links lying in the outside subset ending on 
it. Therefore the operator Qy t . eq.(2.2), which generates gauge transformations at such 
a vertex, maps T-L^ ut to itself. The subset of states in which are gauge-invariant 

with respect to all gauge transformations acting on the outside vertices will be denoted by 
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'H^inv,out ■ Since \iJj) is gauge-invariant, it is easy to see that only states lying in 'h^ inv out 
contribute to the trace in eq.(2.13). This follows from the fact that Qy i has eigenvalues 
±1, and only states lying in T-0^ nv out have eigenvalue +1 with respect to all outside gauge 
transformations. Thus we get 


k 

Pin 


= Tr 


nj k 

ginv ,out 




(2.15) 


The operators p\ n act on 7~Lin- We can define the total electric flux carried inside at a 
boundary vertex, in an analogous fashion, to be the product of the electric flux operators 
for the inside links which end at this vertex. For example, for vertex V\ which has inside 
links 3,4, ending on it the total electric flux operator is 


a 


(3) _( 4 ) 


1 


<7 


1 


(2.16) 


And we can define a vector k, now specifying the total electric flux carried in at each 
boundary vertex. Like 7i 0 ut , the inside Hilbert space, Hj n , can also be written as a direct 
sum, 

U m = ($Uf n (2.17) 

Further, it is easy to see then that since the state \i/j) satisfies Gauss’ law, eq.(2.3), at 
each boundary vertex, the total electric flux going out must equal the total electric flux 
going in at each boundary vertex. This means that the operator pf n will have non-zero 
matrix elements only for states which lie in with k = k. As a result the entanglement 
entropy, eq.(l.l), takes the form, 

See = ~Y^ Tr n k pt, log {pf n ). (2.18) 

k 

As a final step we note that since \ip) is invariant under gauge transformations acting on any 
inside vertex, only states which are gauge-invariant with respect to such an inside gauge 
transformation can contribute in the trace on the RHS above. Denoting by 'H^ inv in the 
subset of all such gauge-invariant states in H^ n , we get, 

See = ~^2t^ pf n log{pf n ). (2.19) 

k 


We are now ready to compare with the results of CHR. Their approach deals with 
only gauge-invariant states and in particular with the algebra of gauge-invariant operators 
acting on these states. Their central point is that this algebra of gauge-invariant operators 
acting on the inside links, Ay, has a non-trivial centre. Once the centre is diagonalised, 
each gauge-invariant subspace corresponding to a fixed set of eigenvalues for the centre does 
factorise, and the entanglement entropy can be defined by tracing over the outside Hilbert 
space. 

Various choices can be made for the centre. One choice, called the electric centre choice, 
corresponds to working with eigenstates of the total electric flux operator entering at each 
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boundary vertex. We see that the different sectors, labelled by the vector, k, in terms of 
which we have expressed our result in eq.(2.19) corresponds to exactly this choice. In fact 
it is easy to see that our result in eq.(2.19) exactly matches the result in eq.(27) of CHR, 
see also eq.(3.12) of [2], for the entanglement entropy. This follows after noting that we can 
define an operator pf n , of unit trace, as, 

p\n = Pkpfn, ( 2 - 20 ) 

with 

Tr^. Pi n = Pk (2-21) 

ginv ,vn 

Eq.(2.19) can then be rewritten in terms of pf n as, 

See = - ^Pklog(pk) + '^2 'PkS{pf n ) (2.22) 

k k 

where S(pf n ). the von Neumann entropy associated with pf n , is 

S(f&) = ~ Tr W k , pf n log(pf n ). (2.23) 

ginv ,m 

Let us end this section with one comment. The different sectors, labelled by different 
values of the electric flux, k, are actually different superselection sectors, as emphasised in 
CHR. Gauge-invariant operators acting on the inside links, or outside links, alone, cannot 
change this electric flux and therefore must preserve the sector to which a state carrying this 
flux belongs. We will return to this point in section 6 when we discuss some measurement 
aspects of the entanglement entropy. 


3 U( 1) Lattice Gauge Theory 

Our discussion for the Z 2 case readily extends to the Z^r and 17(1) cases. Here we will focus 
on the 17(1) case. Once again we work in the Hamiltonian formulation and consider the 
theory on a spatial lattice which for definiteness we take to be a square lattice in 2 spatial 
dimensions. 

Associated with every link on the lattice in this theory is an angular degree of freedom. 
One difference is that the links are directed vectors in this case. The angular degree of 
freedom associated with the link Lij which starts from the vertex i and ends in the vertex 
j is dij E [0, 27t]. With the oppositely oriented link Lji we associate the angle Oji = —Oij. 

Under an infinitesimal gauge transformation at the vertex Vi, parametrised by Cj, the 
angular variables associated with all link variables emanating from V) are shifted by e$, 

Oij —> 0^ + e*. (3-1) 

In the quantum theory associated with each link, L t j, is a Hilbert space, 'Hij- The 
angular variable is promoted to an operator 6^ acting on Hij. There is also an operator 
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[Cij,9jj\ — i 


Cij which is conjugate to Oij , 


(3.2) 


which can be thought of as the angular momentum associated with Oij. It generates shifts 
in the angle 0 tJ : 

( u ‘~'“0j j (: = Oij + e. (3.3) 

The eigenvalues of Cij are quantised since 0ij is compact and takes values over the integers, 
[—oo,.., —1,0,1,oo]. The operator Cij is the electric flux associated with the i — j link. 

An infinitesimal gauge transformation Gvi( e ), parametrised by e, is generated by the 
operator 

g v . = (3.4) 

where the sum is over all links leaving the vertex V). Gauge-invariant states and operators 
meet the conditions, eq.(2.3), and eq.(2.4). 

An underlying Hilbert space CL is obtained by taking the tensor product of the Hilbert 
spaces CLij defined at each link. The space of gauge-invariant states can be naturally 
embedded in CL, CL g i nv C CL. 

For a pure state \C) and a collection of links which we call the inside ones, our proposal 
for the entanglement is completely analogous to the Z 2 case. CL 0 ut,CLi n , are the Hilbert 
spaces given by the tensor product of the outside and inside links. We regard \ijj) £ CL, via 
the natural embedding discussed above, and define pi n as in eq.(2.5). The entanglement 
entropy is then given by eq.(l.l). 

We can show, just as in the Z 2 case, that this entanglement can be written as a sum 
of terms each of which arises from a sector of definite outgoing electric flux. This allows 
us to establish that our definition is the same as the one given in CHR and also [2], We 
adopt the terminology introduced above for the Z 2 case, for inside, outside and boundary 
vertices. The one difference is that we keep track of the orientation of the links, since they 
are directed vectors in the U( 1) case. We take the links which touch a boundary vertex to 
all be oriented to emanate from the vertex. To carry out the sum over CL ou t we then go to 
a basis where the outside links emanating from a boundary vertex are eigenstates of the 
electric flux operators, Cij. The total electric flux leaving a boundary vertex and going out 
is then the sum of the eigenvalues of all the outside links emanating from this vertex. 

The Hilbert space CL ou t splits up into a direct sum of sectors, eq.(2.9), each of which 
is specified by a vector, k, now with integer entries. The i th entry is the total flux leaving 
the boundary vertex Vi. The inside density matrix can then be written as a sum, eq.(2.12), 
with p\ n being given by eq.(2.13). We can also define CL^ inv out C CL^ ut to be the subspace of 
states in flux sector k which are gauge-invariant with respect to all gauge transformations 
on the outside vertices. Then since only states belonging to CL^ inv out contribute to the trace 
in eq.(2.13) we get eq.(2.15). 

For the inside Hilbert space CL ln we now specify the sectors by the vector k so that the 
i th entry is minus of the total electric flux carried by the inside links emanating from the 
Vi boundary vertex. It follows that CLi n can also be written as the direct sum, eq.(2.17). 
Since the state | ip) satisfies Gauss’ law at the boundary vertices, we learn that pf n has 
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non-zero matrix elements with states in T~i\ n only when k = k, leading to eq.(2.18). And 
finally defining T~L^ inv in to be the subspace of states which are gauge-invariant with respect 
to gauge transformations on inside links, we can argue that the entanglement entropy is 
given by eq.(2.19). 

This result is in agreement with the definition given by CHR, see also [2], The different 
sectors which appear in the sum in eq.(2.19) have different values for the electric flux, so 
we see that the result agrees with the electric centre choice of CHR. 


4 The SU(2) Gauge Theory 


The discussion above in turn can be generalised to Non-Abelian gauge theories. As a 
concrete example we take an SU( 2) gauge theory, without matter, and again work in the 
Hamiltonian framework, in 2 spatial dimensions. The degrees of freedom then correspond 
to a 2 X 2 matrix Uij(d) E SU( 2) associated with each link (more abstractly a group 
element of SU( 2)). We see that the matrix is parametrised by three angles 9 which can be 
thought of as the configuration space of a rigid rotor, see [18]. To avoid confusion, let us 
emphasise that the i,j , indices in Uij refer to the link Lij of the spatial lattice with which it 
is associated, and are not the matrix indices of the 2x2 matrix which we have suppressed 
in our condensed notation. The matrix U^- = Uji will be associated with the oppositely 
oriented link denoted by Lji. 

Under an infinitesimal gauge transformation at the vertex Vi the matrix variables cor¬ 
responding to all links emanating from Vi transform as 

Uij^ (l + ie a ^jUij, (4.1) 

where e a = 1,2,3 are infinitesimal parameters, and a a are the Pauli matrices. 

In the quantum theory there is a Hilbert space T~Lij of states associated with each link. 
The matrix variables Uj are promoted to operators, Uij, which act on Uij. Eigenstates of 
Uij will be denoted by | Uij) and provide a basis of Uij. We have 

Uij\Uij) = Uij\Uij) (4-2) 


where Uij is the corresponding matrix in SU{2). 

There are also operators ./“■, a = 1,2,3, which generate SU{2) transformations: 


if a fa ~ _j f b fb 

e vUije le = 


= 1 -ie a 


a 


Uij , 


(4.3) 


where e“ are infinitesimal parameters. On the RHS above these matrices act on the left. 
Acting on the eigenstates, U^) we get, 



(l + ie a 




(4.4) 
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Similarly one can define operators, denoted by J7)“, which result in an infinitesimal group 
element acting on the right, 


e J v\ Uij) = 


Ua ( 1 - ie a< ^— 


(4.5) 


(note the minus sign on the RHS). 

The operators , for each link, L t j. satisfy the standard angular momentum algebra, 

(4.6) 

Similarly, from eq.(4.5) we hnd that the operators J a satisfy the algebra, 

[J a , J b ] =ie abc J c , (4.7) 

see, e..g., [20]. Also, from eq.(4.4), eq.(4.5) it also follows that J a , J h , a, 6, = 1, 2, 3 commute, 

[J a ,J b ]= 0. (4.8) 


Next, consider the vertex V) and orient all the links which touch this vertex to be 
emanating from it. Then a gauge transformation at this vertex is generated by the operator 

g Vi = e i£a ^ J Z (4.9) 

where the sum is over all the 4 links which emanate from V.. Gauge-invariant states \ijj) 
and operators O will then meet the conditions, eq.(2.3) and eq.(2.4). 

It is worth mentioning that the angular momentum operators </?• are not the direct 
analogues of the electric flux operators here. Instead the electric flux carried by the L,j link 
is measured by the difference, 

Qfj = Jij - 3ir ( 4 - 10 ) 

see [18]. This difference arises because gluon degrees carry charge in the Non-Abelian case. 

Consider the directed link Lij which emanates from the vertex Vi and ends on Vj. 
The matrix valued degree of freedom Uij which lives on this link is acted on by both the 
operators J a and J a , which commute with each other, eq.(4.8). The operator J a enters in 
the gauge transformation Qy. eq.(4.9). Similarly the operator J a will enter in the gauge 
transformation at the vertex Vj. This follows after noting that J a and J a generate the left 
and right actions on Uij , eq.(4.4), eq.(4.4) respectively. 

Let us also mention in passing that we can define operators (U(i. t j)) ab . a. b = 1,2,3, 
which acting on the eigenstates | Uij) of the Hilbert space H t] on link L t] give 

(Ug, ij) r b \U i j) = Uft ij) \U ij ) (4.11) 

where is the matrix element of the spin 1 (3-dimensional) representation of SU( 2) 
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corresponding to group element Uij. One can then show that 

J a = -J b U b ^ y (4.12) 

The state on the link Lij can then be specified by the quantum numbers, J 2 = j(j + 
1), J 3 5l 7 3 , for example, since these operators commute with each other. Note that the 
choice of the J 3 ,^ 3 components is purely a matter of convention here; we could have 
equally well chosen other components instead. 

The extended Hilbert space we work with for the entanglement entropy is the tensor 
product of the Hilbert spaces on each of the links, H = (Sl'Hij, Similarly 1-Lout and Hi n 
are the Hilbert spaces obtained by taking the tensor product of the outside and inside link 
Hilbert space. As before, for this case too, we obtain a density matrix, pt n by taking the 
natural embedding of the state |i p) 6 H and then tracing over H ou t , eq.(2.5); and define its 
entanglement entropy to be eq.(l.l). 

It is interesting to try and express our result as a sum over suitably defined sectors 
in this case as well. As before we will use the term boundary vertex to denote a vertex 
that has some inside links and some outside links touching it. At each boundary vertex we 
orient all links to be emanating from it. And define the total outgoing angular momentum, 
J<mt T■ f° be the sum of the angular momenta of all the outside links at this vertex. E.g., 
at vertex V\ shown in Fig 1. Let the two outside links be denoted as 1,2, with angular 
momentum Jf, respectively, then 


= (4-13) 

The states of the outside links can be organised into representations of definite outgoing 
angular momentum, ( J 0 ut,T ) 2 - And the Hilbert space, H 0 uti can be written as a sum over 
sectors of definite values for ( J ou t,T ) 2 - This can be done at each boundary vertex, which in 
general could have 1, 2 or 3, outside links emanating from it, resulting in a total specification 
of the different sectors by giving the quantum number at each boundary vertex for ( J 0 ut,T ) 2 , 
where J° ut T at each vertex is the total angular momentum carried by the outside links. 
The resulting sectors are therefore specified again by a column vector, k, with the entry 
specifying the value of ( J 0 ut,T ) 2 at the boundary vertex, V). 

In a similar way %i n can also be expressed as a direct sum over sectors. We define 
the total ingoing angular momentum, j? T at a vertex to be the total angular momentum 
carried by the ingoing links which emanate from a boundary vertex. E.g., for Vi, with the 
ingoing links labelled as 3,4, we have, 

J 3 + J 4 = Jin } T- (4-14) 

Then the Hilbert space Hi n can be written as a sum over sectors of definite (Jm,r) 2 . And 
Gauss’ law tells us that for the state \i/j) 

( Jout,T ) 2 = («4m,r) 2 (4-15) 

Jout,T W) = -Jin,T IV>),a = 1,2,3. (4.16) 
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We can therefore specify the different sectors which appear in Ki n also by a vector k whose 
entries now stand for the values of ( Ji n ,T) * 2 at the boundary vertices. 

It is worth noting that the difference between the operators J a ,J a , is important in 
ensuring that Ki n can be expressed as a direct sum over different sectors 2 . For example, 
consider the link labelled as “3” in Figure 2, oriented to emanate from vertex V\ and end 
on vertex Vi- The operator J 3 will then enter in the definition of the ingoing angular 
momentum at vertex V\, as discussed above, whereas the operator J£ will enter in the 
definition of the total ingoing angular momentum for vertex V 2 . These commute, eq.(4.8), 
thereby ensuring that Ki n can be expressed as a direct sum over sectors. More generally, 
consider the link Lij which goes from Vi to Vj and the oppositely oriented link Lji. With 
orientation Lij we have operators J?-, 77)“, etc. Then, it follows from the definitions, eq.(4.4), 
and (4.5), that 

= 0 = 1,2,3. (4.17) 

Before proceeding, let us note that a more complete description of the different sectors 
of Ki n is possible, by specifying both ( Ji n ,T ) 2 , and hi addition one component, say Jf n T , 
at every boundary vertex. However, as we will discuss below, this is not necessary because 
fjf n will end up commuting with .J° n Tl and thus will have no non-trivial dependence on 
this extra quantum number. This also makes it superfluous to keep track of additional 
subsectors in K ou t which arise for a fixed value of ( J 0 ut,T ) 2 , and different values of J° ut T . 

The Hilbert spaces, K ou t,'Ki n , can then be written as a direct sum over different sec¬ 
tors of total outgoing and incoming angular momentum at each boundary vertex, eq.(2.9), 
eq.(2.17). The rest of the arguments are essentially identical to the earlier cases. The 
trace over K ou t splits up into different subspaces K^ ut and in each subspace only the gauge- 
invariant states (i.e. states which are gauge-invariant with respect to gauge transformations 
acting on outside vertices) contribute. For a sector specified by the vector ki, say, the re¬ 
sulting operator, [fa , when acting on Ki n only has non-zero matrix elements in the subspace 
77 ^ 1 , dug the gauge-invariance of the state | ip) we started with. Moreover, the operator 
does not have any nonzero matrix elements with states which are not gauge-invariant with 
respect to gauge transformations defined on the inside vertices. From these observations it 
then follows that the resulting entanglement entropy can be written in the form given in 
eq.(2.18) and eq.(2.19). 

To return to a comment made above, it is easy to see, from Gauss’ law, eq.(4.16), and 
the invariance of the trace under a change of basis, that pf n satisfies the relation, 

[Jin,T,pt)= 0, (4.18) 

at every boundary vertex, and for all k. Thus [jf n behaves like a singlet operator under 
JF T , and has no nontrivial dependence on the J^ nT quantum numbers within a sector of 
fixed k. 


2 In general this is true for Hout also, although not for the specific set of inside links considered in Figure 

2 . 
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4.1 Discussion 


We see that our definition yields an expression for the entanglement entropy which is similar 
in form to that obtained in the Abelian cases discussed earlier. 

A few differences are worth taking note of though. The different sectors over which 
the sum must be carried out do not correspond to specifying the total electric flux carried 
away at the boundary vertices, by the outside links, as was mentioned above, since the 
operators J° ut T do not directly measure this flux. Also, we argued above that we do not 
need to keep track of different subsectors specified by the additional quantum number, 
jf nT , say, in each sector with fixed ( Ji n ,T ) 2 , as far as p\ n is concerned, since eq.(4.18) is 
met. However, these subsectors do get important when we discuss more operational ways 
to measure entanglement and compare them with our definition, see section 6. In fact each 
subsector, corresponding to a different value of jf Tl is a different superselection sector by 
itself , which cannot be changed by acting with gauge-invariant operators on the inside links. 
This is because any such operator must commute with ,7" , T . and Gauss’ law, eq.(4.16), 
relates J* ut T to ,Jf n T . 

It is also worth commenting on our result for the Non-Abelian case vis-a-vis the ap¬ 
proach of CHR, [1], which has been discussed above. For the Abelian case the entropy 
obtained with our definition agrees with that obtained in the electric centre choice defini¬ 
tion of CHR. The different sectors, specified by different values for k, which arise in the 
sum in eq.(2.19), correspond to different values for the electric flux operators, which are a 
choice for the centre of the algebra of gauge-invariant operators, Ay. And in this way our 
result directly maps to that obtained in [1] using the electric centre choice. 

In the Non-Abelian case things are similar, as we have seen above. Eigenvalues of the 
operator, ( J 0 ut,T) 2 i or equivalently ( Ji n ,T) 2 , label the different sectors in the expression 
for the entanglement entropy. And this operator is gauge-invariant, although it does not 
directly correspond to the electric flux. In this way our result can be related to a suitable 
choice of a centre in the algebra of gauge-invariant operators, Ay, as suggested by CHR. 

For specifying the different superselection sectors, mentioned above, one can no longer 
work with just gauge-invariant operators. Instead one also needs to specify say the jf n T 
quantum numbers, which are not gauge-invariant. In the extended Hilbert space, 7-/ m , 
which includes non-gauge-invariant states, one can define a bigger algebra containing also 
non-gauge-invariant operators. And the operator jf T , along with (Jj n t) 2 (which is gauge- 
invariant) can be taken to he in the centre of such a bigger algebra. In this way, the 
discussion of superselection sectors could also be perhaps related to the general approach 
of CHR, but now extended to apply to operators in the bigger Hilbert spaces H, 'Kin- 

5 Replica Trick 

We will now show that our formula for the entanglement entropy agrees with a definition 
of the replica trick. Connections to the replica trick are also discussed in [13]. 

The replica trick is based on carrying out a Euclidean path integral in space-time. For 
the (2+l)-dimensional theories we have been studying so far in the Hamiltonian formulation, 
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the path integral is over 3-dimensional Euclidean manifolds. To define the Renyi entropy 
R n = log Trp” n one calculates a path integral over a manifold where the degrees of 
freedom are not single-valued in M 3 but on its ra-fold cover. This n-fold cover is obtained by 
introducing a branch cut at a particular instant of “time,” along the boundary of the spatial 
region of interest, and identifying the values for the variables below the cut in one sheet 
with the values they take above the cut in the adjoining next sheet. A definition for the 
entanglement entropy is then given by suitably continuing the result for the Renyi entropy 
for n — > 1. 

We can apply this procedure readily to the lattice gauge theories we have been studying 
in this paper. The gauge theories have been analysed in the Hamiltonian formulation, living 
on a spatial lattice with continuously varying time. The path integrals which arises for such 
theories most naturally, in the context of our present discussion, are also in the continuum 
in time and on a lattice in the space directions. The degrees of freedom in these path 
integrals are link variables, living on the spatial links, and varying along the time direction. 
E.g., in the U{ 1) case the link variable is an angle, 0 tJ (t) G [0, 2i r], as explained in section 
3, while in the SU(2 ) case it is a 2 x 2 matrix Uij{t) G SU{2), see section 4. 

Now consider a set of links, e.g., the set of of all the links inside and on the boundary 
of the rectangle shown in figure 1. To calculate the entanglement entropy of this set of links 
with the rest, one introduces a cut, say at t = 0, so that in the path integral the values 
these link variables take just above the cut, at t = 0 + , are different from their values just 
below the cut, at t = 0 _ . Then by calculating the path integral on the n-fold cover, and 
taking the appropriate limit, as mentioned above, one obtains the entanglement entropy. 

We will now argue that the result obtained in this way for a gauge theory in fact agrees 
with our definition of the entanglement entropy given in the previous sections. 

5.1 General Background 

To see this, let us first go back to the general case and understand why the path integral on 
the ?i-fold cover is expected to give the Renyi entropy R n . We consider a situation where 
the cut is introduced at t = 0, as in the gauge theory case above, along the boundary of 
the spatial region of interest. The important point is that path integral from (—oo,0 _ ) 
essentially prepares an appropriate state at t = 0 _ , which we denote by \ijj}. Similarly the 
path integral from (0 + ,oo) gives rise to (i/j\. This follows by relating the path integral to 
the evolution induced in Euclidean time due to a Hamiltonian, in the standard fashion. 
With standard boundary conditions, the lowest energy state dominates, resulting in \ijj) 
being the vacuum, 10), of the system. 

For example, take a scalar held theory (which could be on a spatial lattice to make the 
analogy with the gauge theory case more direct). And let us choose a basis of eigenstates 
of the held operator </>(x, t) which we denote as |0(x,f)). Then we get that 

/•0(x,O-) 

(</>(x,f = 0“)|V>) = / [£>0]e (5.1) 

J t= — oo 

where S is the action and we are not being very explicit about the initial conditions at 
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t —>• —cx). A similar expression relates (ip\<p{x, 0 + )) to a path integral from t = 0 + to oo, 

f‘t=00 

(V’I'/HX] t = 0~*~)) = / [ Dcj)\e~ s . (5.2) 

7</>(x,0+) 

The full path integral involves integrating over </>(x, t) while taking the values of the field 
variables at t = 0 to be continuous in the outside region, i.e. 

</>(x, 0~) = 0(x, 0 + ) = 0(x, 0), x E outside, (5-3) 

and taking the values of </>(x, t) in the outside region to be discontinuous, taking values, 
</>(x, 0 _ ) and 0(x,0 + ) below and above the cut. It is then clear from eq.(5.1) and eq.(5.2) 
that this gives rise to the matrix element, 

rt=+oo 

= cr)|p|0(x,f = 0 + )) = / [Dcp]e~ s , (5.4) 

J t= — OO 

with x E inside. In eq.(5.4) p is the density matrix of the inside region, and the path 
integral is carried out in the presence of the cut. 

Once this identification between the path integral on one sheet and the matrix element 
of the density matrix is established, it follows immediately that the Renyi entropy R n is 
given by the path integral over the n-fold cover. And the entanglement entropy can then 
be calculated from the formula, 

*5 replica -[( nd n - 1) logTrp n ]| n= i. (5.5) 

assuming the continuation to n —> 1 is well defined. A minor point is that the path integral 
does not usually give rise to a normalised state \ip) of unit norm, but eq.(5.5) takes this 
into account to give the von Neumann entropy of the correctly normalised density matrix 
with unit trace, see, e.g., [9]. 

5.2 Gauge Theories 

We can now return to the case of gauge theories. For this case a basis which is the analogue 
of | cj)(x,t)} above can be introduced as follows. We start with eigenstates of the operator 
Uij at each link. The operator JJij for the SU( 2) theory was introduced in section 4 and 
its eigenstates were denoted by | Uij), with the eigenvalues U t j being elements of SU( 2). 
Similarly, an operator Uij can be defined in the 17(1) case with its eigenvalues Uij being 
elements of U{ 1). A full basis can then be obtained by taking the tensor product of all such 
basis elements on each link. An element of this basis is given by | Uij). This basis 

is the analogue of |0(x, t)). A basis for the inside links is obtained by taking the tensor 
product only over the basis elements for the inside links, which we schematically write as 
®(v)e*n | Uij). Similarly a basis for the outside is given by I Uij). 

It is then easy to see that the discussion above applies also for these gauge theories. The 
path integral from (—oo,0 _ ), and (oo,0“ l "), gives rise to the state |'i/>) and (ip\ respectively. 
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With eq.(5.1) taking the form, 


r Uij (o ) 

(Uiji0~M = / [DUij] e~ s (5.6) 

J t= — oo 

and similarly for eq.(5.2). For the 51/(2) case the functional integral for Uij on the RHS of 
eq.(5.6) involves the standard Haar measure on SU{ 2). 

The full path integral is obtained by also integrating over the outside link variables at 
t = 0, we denote these by Uij(t = 0) with (ij) € out . These variables are continuous across 
t = 0. This gives rise to the density matrix p acting on the inside links. As a result, the 
path integral with U^- )&nslde taking values, at t = 0" and U( ij ) &nside (0 + ), 

at t = 0 + respectively, leads to the matrix element 

(U(ij)&n{® + )\p\U(ij)£in{ 0 )) = j [DUij\e s . (5-7) 

It then follows that the n-fold cover computes the Renyi entropy Trp n and thus this defini¬ 
tion of the replica trick can be interpreted as the von Neumann entropy for p. The action 
S is related to the Hamiltonian in the standard way by replacing the momenta by time 
derivatives of Uij. 

Note that the path integral on the RHS preserves gauge-invariance with respect to 
spatial gauge transformations, i.e. gauge transformations where the parameters vary only 
in space and not time. This follows from the invariance of S and the invariance of the Haar 
measure under these gauge transformations. If the initial state at t = — oo is also gauge- 
invariant, or in effect gauge-invariant as would happen if the vacuum dominates, then \ijj) 
at t = 0“ would also be gauge-invariant and one can show from eq.(5.6) that 

(UiM) = m\1>) (5.8) 

where Uij, Uf- are related by the gauge transformation g. 

Usually the path integral formulation for lattice gauge theories is presented in a some¬ 
what different manner, corresponding to the Euclidean formulation, with both time and 
the spatial coordinates being on a lattice. Starting from this formulation one can obtain 
the above type of path integral representation by taking the continuum limit along the t 
direction and fixing Aq = 0 gauge. These steps are identical to those taken in obtaining 
the Hamiltonian formulation. The resulting theory then only has invariance with respect to 
spatial gauge transformations. Note that no further gauge fixing is needed to make the path 
integral in eq.(5.7) well defined, since the integrals involve the Haar measure on compact 
groups and are finite. 

We are now ready to show that the replica trick definition given above agrees with our 
definition of the entanglement given in the previous sections. In fact to do so we simply have 
to note that the discussion above has, in effect, been carried out with \ij)) being embedded 
in the extended Hilbert space, H. This extended Hilbert space, we remind the reader, 
includes in particular non-gauge-invariant states. This is because the basis ®(ij\ \ Uij) is 
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in fact a basis of 'H and similarly <S)(vj')em Wij) and <S>(ij)£out Wij) are a basis for Bin 
and 1-Lout respectively. Thus the density matrix defined by the path integral is actually 
Pin = Tr% otit (|'0)( , 0|), eq.(2.5), and the Renyi entropies and von Neumann entropy obtained 
from the replica trick agree with those obtained by taking the trace of suitable powers of 
Pin in 'Hi n . This shows that the answer one gets from the replica trick defined above agrees 
with our definition in eq.(l.l) for the entanglement entropy. 

We have considered the SU{ 2) and f7(l) cases above in 3 space-time dimensions, for 
definiteness, but it is clear that the argument applies more generally, in higher dimensions, 
and also for other Non-Abelian gauge groups, and Abelian discrete groups. 

Let us end this section with one comment. On the lattice other definitions of the 
replica trick are possible, different from the one we have used above. One expects that in 
the continuum limit all these definitions agree unto counter terms, including some localised 
on the boundary of the region of interest. 3 

6 More on The Entanglement 
6.1 Measuring the Entanglement 

ft is well known in the literature of quantum computation that the entanglement in a 
bipartite system can be quantified operationally by comparing it to the entanglement of a 
reference system, usually taken to be a set of N Bell pairs. 4 This comparison can be done 
in two ways, see, e.g., [21], [22], For example, in the process of entanglement distillation, 
one starts with the entangled parts A, B, of a bipartite pure system and an extra set of 
unentangled qubits. The two subsystems A, B , are assumed to be separated in space. 
We then ask, in an asymptotic sense with a large number of copies of the initial system, 
how many pairs of entangled Bell pairs can be produced by transferring the entanglement 
between A, B to the qubits, using only local operations and classical communication. For 
a system with local degrees of freedom, e.g., a spin system, it is known that this number N 
is given by 

N = —Tr(plogp) (6.1) 

where p is the density matrix of A or of B. Alternatively, in entanglement dilution, we take 
the two parts A, B to be nearby, and a reservoir of Bell pairs such that one partner in each 
pair is near the system AB, while the other is far away. We then ask how many Bell pairs 
would be needed for teleporting one of the two parts, say A, to the distant location. The 
answer, for the number of such Bell pairs, N, in the asymptotic sense mentioned above, for 
a system with local degrees of freedom is again given by eq.(6.1). 

We would now like to examine whether for gauge theories under consideration here the 
result of such a comparison with Bell pairs also agrees with our definition of the entangle¬ 
ment entropy given in eq.(l.l). In fact, this turns out not to be true. The entanglement 
entropy which we have defined is not equal to the number of Bell pairs involved in entan- 

3 We thank Gautam Mandal and Shiraz Minwalla for a discussion on this point. 

4 We thank John Preskill for insightful comments which led us to examine issues discussed in this sub¬ 
section. 
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glement distillation, or dilution, in general. The two parts A, B of the system, are the set 
of inside and outside links, as per the definitions used above. The essential obstacle is that 
the local operations which one can use correspond to gauge-invariant observables acting on 
the inside or the outside links alone. We saw above that our result for the entanglement 
entropy can be expressed as a sum over different superselection sectors. Gauge-invariant 
local operators cannot connect these different sectors, and acting with only these operators 
is then not enough to make the full entropy available for a comparison with Bell pairs. 5 

In the Abelian case the different superselection sectors correspond to different values 
of electric flux, and this observation about the full entropy not being available for the 
comparison was also made by CHR. Working within a flux sector labelled by k the entropy 
can be converted into entangled Bell pairs, or alternatively used in teleportation. But the 
rest of the information in See , eq.(2.22), having to do with the probabilities p ^ is not 
directly related to the standard measures of entanglement. Note that a dependence on p ^ 
enters both in the first term and also the second term in eq.(2.22) through the coefficients, 
Pki so neither of these terms or their sum is directly the result of such a comparison. 

In the Non-Abelian case, as was mentioned in section 4, the different superselection 
sectors do not correspond just to different values of k, and therefore to different values of 
( Jin,T ) 2 - Rather, the different superselection sectors are given by subsectors within each 
sector of fixed k. These subsectors correspond to different values of a component of the in¬ 
coming total angular momentum, say Jf n T , for a fixed value of (Jm,T) 2 7 at each boundary 
vertex. It would be worth understanding this difference between the Abelian and Non- 
Abelian cases further. 

It is also worth mentioning that even though our answer for the entanglement, eq.(l.l), 
does not directly correspond to the number of Bell pairs which are obtained in entanglement 
distillation or dilution, the density matrix pi n we obtain can be completely measured by 
local experiments performed on the inside links alone, as long as one has an ensemble of 
many copies of the inside system available to us. For example, in the Abelian case, one can 
measure the electric flux on each of the links, and also transitions induced between different 
values of these electric fluxes by inserting Wilson loops lying in the inside, and these are 
enough to determine pi n completely. 

The situation for the Non-Abelian case is similar. We saw in eq.(4.18) that p\ n is a 
singlet under the total incoming angular momentum, ,J" n T , at each boundary vertex. For 
determining p\ n in each sector, from which the full density matrix can be obtained, it is 
therefore enough to obtain the values of Tr (pf n O) for a complete set of Hermitian operators 
which are singlets under J? T . In fact gauge-invariant observables acting on the inside links 
provide such a complete set, and therefore allow us to determine p m fully. 

’In gauge theories it is still true that for a large number of copies, n, of the system of interest one can 
work in effect in a subspace of reduced dimensionality, 2 P , where p = uSee, and See is the full entanglement 
entropy given by eq.(l.l). This smaller dimension subspace also admits a tensor product decomposition 
between A and B. But the process of distillation requires a further unitary transformation which transfers 
the entanglement between A,B, in this reduced subspace to a set of Bell pairs. For gauge theories this 
unitary transformation in general cannot be carried out using gauge-invariant operators alone, since these 
operators cannot connect different superselection sectors. As a result, the full entanglement entropy we 
have defined, eq.(l.l), is not available for distillation. 
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Figure 3. The red curve is the boundary of the spatial region of interest. Several links intersect 
it, including link A which emanates from vertex V 4 and and links F,G, which emanate from vertex 
Ve- 


To see why gauge-invariant operators acting on the inside are singlets we note that a 
gauge-invariant operator acting on the inside must commute with the total outgoing angular 
momentum at each boundary vertex, J° ut T , by definition because the operators only act 
on the inside links. And they must also commute with the total angular momentum, 
Jj, = Jf n T + since these are the generators of gauge transformations. Thus they 

must commute with ,J" n T . which means they are singlets with respect to these operators. 
A little more thought shows that they are also a complete set of such singlet operators. 

6.2 Entanglement For Spatial Regions 

In the introduction, we began by considering a region of space and asking about its entan¬ 
glement, but then chose to focus on the entanglement for a set of links instead. Now that 
we have understood the situation for a collection of links, to some extent, it is worth going 
back and asking about the entanglement for a region of space. Suppose for concreteness 
that we are interested in the region of space enclosed by the curve in Fig 3, which has one 
connected component. There are links leaving this region, which intersect the curve; these 
are among the dashed links shown in Fig 3. To use our definition for the entanglement 
entropy in such a situation we need to decide whether to include the intersecting links in 
the set of interest (the inside as per our definition above). 

To be specific let us take the case of an Abelian gauge theory. We will see below that 
for some links which intersect this curve the answer for the entanglement of the inside set of 
links does not change whether or not they are included. But for other links it can change. 
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Figure 4. The extended lattice construction. Blue dots represent additional vertices which have 
been introduced at the points where links of the original lattice intersect the boundary shown in 
red. Links touching the blue dots are the new links which have been introduced. 


First, consider link A, as shown in Fig 3. This is the only link which connects the 
boundary vertex Va to the outside. Gauss’ law at vertex Va fixes the value of the electric 
flux at this link in terms of the electric fluxes carried by links B, C , D which are inside. 
Similarly Gauss’ law on vertex Va 1 which is outside and the other end-point of link A fixes 
the value of the flux on link A in terms of the flux at the three other links, B’,C’,D’, which 
are outside and also end on Va>- This means that the electric flux on A is completely 
correlated with degrees of freedom inside and outside. It then follows simply that including 
this link variable in the inside or the outside does not change the entanglement entropy 
eq.(l.l). Next, consider the vertex Ve which is connected to the outside by two links F,G, 
and to the other inside vertices by H,J. It is now clear that Gauss’ law at Ve only relates 
the total electric flux on F,G, with that on H,J, and does not fix the individual values of 
the two electric fluxes on F,G. This means there are degrees of freedom on the F,G, links 
which are not fully correlated with the inside and outside (in the Z 2 case this is one qubit’s 
worth). As a result, the entanglement we obtain by including or not including links F, G, 
will be different in general. 

The example above shows that ambiguities arise when we attempt to extend our def¬ 
inition of the entanglement entropy to apply to a spatial region instead of a collection of 
links. 

6.3 The Extended Lattice Construction 

Another definition for the entanglement entropy both in Abelian and Non-Abelian gauge 
theories has been proposed in the literature, see, [11], [12], [13], [14], [17]. Following [1], 
which has a helpful summary that we will use here, we will refer to this as the extended 
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lattice construction (ELC). This definition applies to a region of space whose boundary 
does not go through any vertex of the lattice. In figure 4 we show such a region which 
lies in the interior of the boundary red curve. The boundary is intersected by links which 
connect the inside to the outside. The idea behind the definition is to extend the lattice 
by introducing extra links and vertices such that the enlarged Hilbert space associated 
with this extended lattice then admits a tensor product decomposition. This is done by 
introducing extra vertices at points where links intersect the boundary, and splitting each 
link which intersects the boundary into two by allowing it to terminate at the extra vertex 
that has been introduced. For example, in figure 4 consider the link which connects Va and 
Va'i f° r brevity we refer to it as the A — A! link below. We introduce a new vertex Va" 
where this link intersect the boundary and then split the A — A' link into two, giving rise 
to the A — A" , and A" — A links. The extra vertices introduced in this way are shown in 
blue in figure 4. 

In the ELC definition the Hilbert space associated with the extended lattice, , con¬ 
sists of states which satisfy the Gauss’ law constraint at the original vertices, shown in black 
in Fig 4, but which need not satisfy this constraint on the extra vertices which have been 
introduced, shown in blue. As a result it turns out that %' now admits a tensor product 
decomposition. 

To understand this better, let us introduce the terminology of “inside links” to refer to 
links which either he securely inside the boundary or to links which connect vertices lying 
inside with the extra vertices which he on the boundary. Outside links will then refer to all 
remaining links. For example, in figure 4 the inside links include the link A — A" , but not 
the link A" — A' which lies in the outside. It is then easy to see that since the states in %' 
do not necessarily meet the Gauss’ law constraint on the extra vertices we have the tensor 
product decomposition 

W = H' in ®H' out . ( 6 . 2 ) 

Here 'H , in is the Hilbert space of states defined on inside links which are gauge-invariant 
with respect to inside vertices but not the extra vertices lying on the boundary curve. And 
similarly 'H , out is the Hilbert space of states on outside links which are gauge-invariant with 
respect to the outside vertices but not the vertices on the boundary curve. 

A state | ip) which lies in T-L g i nv - the Hilbert space of gauge-invariant states defined on 
the original lattice - is then embedded as a state |i p 1 ) E 'H! in a well defined way, given 
by eq.(6.3). Although 'H' includes in general states which are not gauge-invariant with 
respect to gauge transformations at the extra vertices, as mentioned above, all states, \ip') 
which arise in 'H' due to embedding in this manner are actually gauge-invariant. As an 
example consider the links A — A ", A" — A', in Fig 4, and let the link variables on them 
take values Ua~A” > Ua v -A’i respectively. Then the state \tp') E W is defined by its wave 
function satisfying the relation 

(Ua-a", Ua"-a') = Wa-a" • Ua"-a')i (6-3) 

where denotes matrix product. Here we are suppressing the dependence of ip on the 
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other link variables. Under a gauge transformation parametrised by group element g at the 
extra vertex A", 

Ua-a" -> Ua~a" ■ g~ x , Ua"~a' —t g ■Ua"~a i - (6-4) 

As a result we see that \ijj') is gauge-invariant. 6 One minor comment is that | ijj') as defined 
in eq.(6.3) does not necessarily have unit norm. This can be taken care of by introducing 
a normalisation constant on the RHS of eq.(6.3), or in the case of continuous groups, like 
the t/(l) or SU(N ) examples, by defining the measure for summing over the group to be 
of unit norm (as is the case with the Haar measure). 

Since Ti 1 admits a tensor product decomposition one then defines the entanglement 
entropy of the inside links by tracing over T~L' out to obtain a density matrix, 

p = Tr K J^)(^\, (6.5) 

and then equating the von Neumann entropy of p, 

S = -Trw in P log(p), (6.6) 


with the entanglement entropy. 

We will now argue that both for Abelian and Non Abelian cases the ELC definition 
above agrees with our definition, eq.(l.l) when applied to the extended lattice. We have 
also argued that in the Abelian case our definition is equivalent to the electric centre choice 
of CHR. It will then follow that for the Abelian case the ELC definition is equivalent to 
the electric centre choice as well, when applied to the extended lattice, as was noted in [1] 
already. 

The equivalence with our definition is in fact straightforward. Consider applying the 
definition we have proposed, eq.(l.l), to the extended lattice described above, with the 
choice of the inside and outside subsets of links, for purposes of applying our definition, 
being the same as what was given above in the ELC case. And with the gauge-invariant 
state with which our analysis begins being the state | if}') on the extended lattice. 1 It is then 
easy to see that in fact 'H ! out , 'H l m , introduced above can be expressed as the direct sums, 


n_/' _ CTA -l j 

n out — n 


k 

ginv^outi 


(6.7) 


and, 


qj> = CTA -i/k 

rL in VL7 rL ginv, 

k 


( 6 . 8 ) 


where the subspaces 'H^ inv out . 'H^ invin were introduced in section 2. Since pi n is given by 
eq.(2.12), with p k being defined in eq.(2.15), we see from eq. (6.7) that pi n is in fact the 
same as p given in eq.(6.5). Similarly the trace in eq.(2.19) can be seen to be the same as 
the trace over thereby leading to the entanglement entropy being the same. 


f, This discussion applies also to the Non-Abelian case, see, [14], for which the order in which g acts in 
eq.(6.4) is important. 

'This means that for applying the equations in section 2, e.g., eq. (2.11), we replace \ip) by | ip'). 
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We end with one comment. The ELC definition pertains to the extended lattice and 
does not directly yield the answer for the entanglement entropy of a set of links in the 
original unextended lattice. For example, consider the set of links shown in bold face which 
lie inside the red curve in figure 4. It is easy to see that the entanglement entropy which 
arises using our definition for this set of links is different from the entanglement entropy for 
the extended lattice associated to it, in general. This difference arises in the Abelian case 
itself, and is also present in the non Abelian case. 

Let us begin with the case of an Abelian theory, and let us consider the situation for 
the A — A 7 link in the original lattice. As argued in the previous subsection, in the original 
lattice, the electric flux carried by link A — A' is correlated completely with the total electric 
flux on the inside links which terminate at Va and also the outside links terminating at Va> ■ 
In the extended lattice, given the mapping proposed in the ELC construction for a state 
in the extended Hilbert space , one can show that the value of the electric flux carried by 
the A — A” and A ” — A links are the same, and equal to the value of the electric flux on 
the A — A' link in the original lattice. It therefore follows that the two new links, A — A” 
and A” — A', of the extended lattice also take values completely correlated with the total 
electric flux carried by the inside links at the Va vertex and the other outside links at the 
Va” vertex respectively. It is then easy to see that adding these two links, A — A”, A” — A', 
to the inside and outside sets does not change the entanglement. However, if we instead 
consider vertex Ve, also discussed in the previous subsection, the conclusion is different. 
Now in the original lattice only the total electric flux on the outside links F, G is fixed by 
Gauss’ law applied to the vertex Ve- The individual fluxes on the links F, G are not fixed. 
As a result, one finds that when both of these link are split into two, and the procedure in 
the ELC is carried out, the entanglement we get will in general be different from the one 
calculated by using our definition (or equivalently the electric centre choice of CHR) on the 
original lattice. 

In the non-Abelian case also there is a difference between the ELC definition and the 
one we propose. In fact the difference is if anything even more pronounced and arises 
even for a link of the A — A' typed discussed above, due to the fact that more quantum 
numbers are needed to specify the state of a link completely in this case. In the SU{2) 
gauge theory for example, the state of the A — A' link is specified by 3 quantum numbers, 
(J) 2 ,J 3 ,J 3 , as discussed in section 4. Gauss’ law at Va fixes ( J ) 2 , J 3 , and at vertex Va>, 
(J) 2 = {J) 2 ,J 3 , in terms of the total angular momentum of the inside links and the other 
outside links emanating from Va, Va 1 respectively. When the link is divided into two, as per 
the definition of the ELC construction, the (J) 2 , J 3 numbers of the A —A” link are fixed by 
the total angular momentum of the inside links at Va and the (J) 2 ,J 3 quantum numbers 
of the A” — A' link are fixed by the total angular momentum of the other outside links 
at the vertex Va>- However Gauss’ law does not automatically determine the J 3 quantum 
number of A —A”, and the J 3 quantum number of A” — A 7 in terms of the angular momenta 
of the inside links on Va and the outside links on Va> respectively. Gauss’ law at the new 
vertex A” does tells us that these two additional quantum numbers must take equal values 
though and they are therefore fully correlated. It is then easy to see that including A — A” 
in the inside and A” — A 7 in the outside will in general gives a result which is different from 
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what is obtained using our definition for the original set of links on the unchanged lattice. 


7 Conclusion 

In this paper we have proposed a definition for the entanglement entropy of a gauge theory. 
Our definition applies to both Abelian and Non-Abelian gauge theories on a spatial lattice. 
In the Abelian case without matter, we show that the definition is equivalent to the electric 
centre choice of CHR, [1]. In the Non-Abelian case, we find that the situation is similar and 
the definition is tied to diagonalising operators lying in the centre of the algebra of gauge- 
invariant operators, Ay, defined on the links of interest. We also show that in general, 
both for the Abelian and Non-Abelian cases, our definition agrees with the result obtained 
by using a definition of the replica trick. This is a good argument for taking the definition 
given here seriously. 

Our definition applies to any subset of links on the spatial lattice and results in standard 
desirable properties, like strong subadditivity. Interestingly, we find that the definition we 
propose, does not agree with standard ways of quantifying entanglement discussed in the 
quantum information theory literature. These ways are based on using a collection of 
Bell pairs as a standard measure of entanglement. An entanglement entropy is assigned 
between two parts of a system (the inside and outside links) by comparing them with a 
collection of Bell pairs, using the process of entanglement distillation or dilution which 
involve local operations and classical communication. The difference between our definition 
and such a more operational definition has to do with the presence of superselection sectors. 
These superselection sectors cannot be changed by local gauge-invariant operations acting 
on only one of the two parts, and this prevents the whole entanglement from being made 
available for the process of comparison. For Abelian theories the different superselection 
sectors correspond to different values of the electric flux, as noted by CHR. In the Non- 
Abelian case this is not true any more. In fact the situation is even more interesting since 
the superselection sectors are labelled by quantum numbers some of which are not gauge- 
invariant. This non-gauge-invariance drops out of the final answer for the entanglement 
which is obtained by summing over all sectors. 

We find this mismatch between our definition and the more operational ways to quan¬ 
tify entanglement quite interesting. At root, the mismatch arises because gauge theories 
have non-local gauge-invariant operators, like electric and magnetic flux loops, and corre¬ 
spondingly non-local excitations. It might therefore be a general feature of attempts to 
define the entanglement entropy in gauge theories. 

There are many other ways in which the work reported here can be extended. This 
includes analysing the continuum limit, other measures of entanglement, like mutual infor¬ 
mation and relative entropy, and a study of gauge theories with matter. There is a good 
argument now showing that the Ryu-Takayanagi (RT) proposal in AdS gravity theories, 
[23], [24], follows from the replica trick in the dual field theory, [25]. The definition we give, 
by embedding the state in an enlarged Hilbert space, H, can be applied to matter theories 
as well. And it seems quite reasonable that the equivalence between our definition and the 
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replica trick would continue to hold in these cases too. In this way one can hope to firmly 
establish the equivalence between our definition and the RT proposal in AdS gravity. 

Let us end by mentioning some more topics for further study. It would be worth think¬ 
ing about the entanglement entropy for a spatial region, as opposed to a set of links, further, 
since preliminary attempts on the lattice at least lead to ambiguities, as was discussed above. 
Also, the behaviour of the entanglement entropy under a duality transformation, which ex¬ 
changes a gauge theory with a spin system in 3 dimensions, for example, is worth elucidating 
further. Finally, the connection to the extended lattice construction, as discussed in 6.3, 
too bears further thought. 
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